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In  this  paper,  we  consider  the  linear  quadratic  optimal  control 
problem  on  infinite  time  interval  for  linear  time-invariant  systems  defined  on 
Hilbert  spaces.  The  optimal  control  is  given  by  a  feedback  form  in  terms  of 
solution  n  to  the  associated  algebraic  Riccati  equation  (ARE).  A  Ritz  type 
approximation  is  used  to  obtain  a  sequence  nN  of  finite  dimensional 
approximations  of  the  solution  to  ARE.  A  sufficient  condition  that  shows  nN 
converges  strongly  to  n  is  obtained.  Under  this  condition,  we  derive  a 
formula  which  can  be  used  to  obtain  a  rate  of  convergence  of  nN  to  n. 
We  demonstrate  and  apply  the  results  for  the  Galerkin  approximation  for 
parabolic  systems  and  the  averaging  approximation  for  hereditary  differential 
systems. 
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Assume  Z,  U  and  Y  are  Hilbert  spaces.  Consider  the  evolution 
equation  on  Z 

(1.1)  z(t)  =  A  z(t)  +  B  u(t)  ,  z(0)  =  z0  €  Z 

where  u(t)  is  a  U-valued  control  function,  A  is  the  infinitesimal  generator  of 
strongly  continuous  semigroup  S(t)  on  Z,  and  B  e  £(U,Z).  The  Y-valued 
observation  function  y  is  given  by 

(1.2)  y(t)  =  C  z(t)  ,  t  >  0  . 

We  assume  that  C  €  £(Z,Y).  We  interpret  the  equation  (1.1)  in  the  mild  sense; 
the  solution  of  (1.1)  is  given  by 

(1.3)  z(t)  =  S(t)z„  +  I*  S(t-s)R  u(s)ds  . 

■*0 

Consider  the  minimization  problem;  minimize  the  cost  functional 

(1.4)  J(u,z0)  -  [  (||y(t)||a  +  ||u(t)  ||a)dt 

Jo 

subject  to  (1.3).  Then  the  following  result  is  well-known  [10],[1I]: 

Theorem  1.1  Assume  (A.B)  is  stabilizable  and  (A,C)  is  detectable.  Then 

there  exists  a  unique  nonnegative  self-adjoint  solution  n  to  the  algebraic 
Riccati  equation  in  Z: 

(1.5)  ( ATI  +  IL4  -  rWB'Tl  +  C*C) z  =  0  for  all  z  €  dom(/l)  . 
and  the  optimal  solution  u°  to  (1.4)  is  given  by 

u°(t)  =  -R*n  T(t)z0 

where  T(t)  is  the  strongly  continuous  semigroup  generated  by  A  -  BB* n 
and  it  is  uniformly  exponentially  stable. 
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Herc  wc  have 
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Definition  1.2  (1)  (A,B)  is  stabilizable  if  there  exists  an  operator 

K  €  £(Z,U)  such  that  A-BK  generates  a  uniformly  exponentially  stable 
semigroup  on  Z. 

(2)  (A.C)  is  detectable  if  there  exists  an  operator.  G  €  £(Y,Z)  such  that  A  -  GC 
generates  a  uniformly  exponentially  stable  semigroup. 


The  purpose  here  is  to  construct  a  finite  dimensional  approximation 
of  the  optimal  feedback  gain  operator  B*Tl  Let  us  consider  a  sequence  of 
approximating  problems  (ZN,/tN,5N,CN);  let  ZN  be  a  sequence  of  finite 
dimensional  subspaces  of  Z  and  PN  be  the  orthogonal  projection  of  Z  onto 
ZN.  Assume  /1N  :  ZN  -*ZN  ,  SN  :  U  -  ZN  and  CN  :  ZN  -»  Y  are  continuous. 
Then  consider  the  Nth  approximating  problem  of  (1.4) 


(1.6)  minimize  J 


N(u,z0)  *  jo  (||CNzN(t)||J  +  ||u(t)||2)dt 


9ft 

$ 


subject  to 


(1.7)  zN(t)  =  SN(t)PNz0  +  [  SN(t-s)ffNu(s)ds 

*  n 


where  SN(t)  =  cA  t  0.  Then  the  optimal  control  uN  of  (1.6)  is  given 


uN(t)  =  .flNVeMN-*N*NnN)tPNZo  ,  t  ?  o 


where  nN  :  -  ZN  is  self-adjoint  and  satisfies  the  Nth  approximating 

algebraic  Riccati  equation  in  ZN; 


(1.8)  a"'  V  +  nN^N  -  nNBNBN  nN  +  cN  cN  =  o 


Here,  BN  n  ,  N  >  1  yields  a  sequence  of  finite  dimensional  approximations  of 


the  optimal  feedback  gain  [3]. 

In  this  paper  we  first  obtain  a  condition  on  (ZN^4N,BN,CN)  for  which 
(1.8)  admits  a  unique  nonnegative  solution  nN,  and  nNPN  converges  strongly  to  n 
in  §2.  Such  a  condition  has  been  discussed  in  [2],  [3]  but  the  condition  in  this 
paper  improves  those  in  [2],  [3],  i.e.,  we  introduce  the  uniform  detectability 
condition  (see,  (H3)  in  §2,  for  the  definition)  which  is  additional  to  those 
considered  in  [2],  and  using  this  condition,  we  are  able  to  show  that  there 
exists  an  integer  N0  such  that  for  N  N0 

jje(/lN-fiNSN  nN)tpN||  4  m  e‘ut  ,  t  >  0 

for  positive  constants  M  *  1  and  u  (independent  of  N  >  N0).  This  assertion 
is  a  part  of  assumptions  in  [2,  Theorem  2.2],  The  uniform  detectability 
condition  is  satisfied  if  C*C  is  coercive,  which  is  assumed  in  the  discussions  in 
[2,p.  693].  Thus,  the  uniform  detectability  condition  can  be  regarded  as  a 
relaxation  of  the  coercivity  assumption  mentioned  above.  Next,  under  the 
condition  in  §2  we  derive  a  formula  which  provides  a  rate  of  convergence  of 
nN  to  n  and  apply  the  formula  for  specific  examples. 
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We  assume  the  following.  Let  SN(t)  =  e^Nt,  t  >  0 


(HI)  For  each  z  e  Z,  we  have 

(i)  SN(t)PNz  -  S(t)z  ,  and 

(ii)  SN(t)*PNz  -  S*(t)z  , 

where  the  convergences  are  uniform  in  t  on  bounded  subsets  of  [O,00). 

(H2)  (i)  For  each  u  «  U,  SNu  -*  fiu  and  for  each  z  e  Z 

Bn*  PNz  -  B*  z. 

(ii)  For  each  z  e  Z,  CNPNz  •*  Cz  and  for  each  y  e  Y 
CN*y  -  C* y. 

(H3)  (i)  The  family  of  the  pairs  C4N,2?N)  is  uniformly  stabilizablc:  i.c. 

there  exists  a  sequence  of  operators  KN  e  £(ZN,U)  such  that 
sup  |{A'N||  <  •  and 

||e(/lN-BN/:N)tpN|j  <  M1e*u)it  ,  t  >  0 
for  some  positive  constants  Mx  >  1  and  u)x  . 

(ii)  The  family  of  the  pairs  (/tN,CN)  is  uniformly  detectable;  i.c. 
there  exists  a  sequence  of  operators  GN  €  *(Y,ZN)  such  that 
supj|GN||  <  •  and 

/  jN  rrN/"N\t  N  -u,t 

||e^  G  c  >lPN||  (  M2e  2  ,  t  >  0 
for  some  positive  constants  M2  )  1  and  u>2. 

Remark  (1)  Suppose  =  PN£  and  CN  =  CPN.  Then  (H2)  holds 

since  it  follows  from  (HI)  that  PNz  -»  z  for  all  z  e  Z. 

(2)  The  assumption  (H3)  is  closely  related  to  the  preservation  of 
exponential  stability  under  approximation  in  [3, Conjecture  7.1]  and  it  is  shown 
in  [2]  that  (H3)  (i)  ((POES)  in  (2])  is  satisfied  for  parabolic  systems  using  the 


Galerkin  approximation. 


I  H* M 


i 


(3)  A  natural  way  to  argue  (H3)  is  to  take  /fN  =  KPN  and  GN  =  PNC7  for 


some  K  e  £(Z,U)  and  G  t  £( Y,Z)  such  that  A  -  BK  and  A  -  GC  generate 


uniformly  exponentially  stable  semigroups  on  Z. 


Theorem  2,1  Suppose  (H1)-(H3)  are  satisfied.  Then  for  each  N,  (1.8) 


admits  a  unique  nonnegative  solution  nN,  supl|nN||  <  ®,  and  there  exist 


positive  constants  Ms  *  1  and  u3  (independent  of  n)  such  that 


(a N  RNRN*nN\*  x,  -u,t 

||e^  -B  B  n  )tpN  ||  i  MsC  3  ^  t  J  o 


Proof:  The  proof  is  based  on  the  arguments  in  [11],  The  existence  and 


uniqueness  of  solutions  to  (1.8)  follow  from  Theorem  1.1.  Since 


<nNPNz,z>  =  min  JN(u,z),  (H3)  (i)  implies  that 


<nNPNz,z>  <  Jn(-Knzn(-);z) 


[  (||CNe(/,N‘fiN*N)tPNz||*  +  ||KNe^N‘^N^N)tPNz  ||2)dt 

*  n 


<  B||z||J  for  some  positive  constant  B 


Since  nN  is  self-adjoint,  nonnegative  definite,  this  implies  that  ||nN||  s  B. 


By  the  variation  of  constants  formula 


(2.1)  e(<4N-#NflN*nN)t  .  TN(t)  +  r‘  TN(t-s)(CNCN-5N5N*nV'4N’BNsN*nN)s  ds 

*  n 


where  TN(t)  =  e^  ^  t  >  0.  Here,  from  (1.8) 


(,4N-flNBNV)*nN  +  nN(/tN-BNBN*nN)  +  nNsNsNV  +  cN*cN  =  0  , 


so  that  if  zN(t)  =  e^-^^VjipHz,  t  >  0,  then 


I 


k 


I 


& 

*> 

v 

‘!i 


i 


ij»; 

!*i< 


ft! 

»v 


B 


y 

V,V, 


Thus,  for  all  t  >  0 


(2.2)  <nNzN(t),zN(t)  >  +  f(||BNVzN(t)||2  +  ||CNzN(t)||2)dt 

«  <nNPNz,z>  S  B||zj|2 

Now,  from  (2.1),  we  have  for  all  t  *  0 

fl  M2  2M2  ft 

J  llzN(s>l|2ds  <  -JT  !lzll2  +  ~2  (»cNil2  +  llfiNll2)  J  nNzN(s>H2 

*  0  2  Uj  J0 

+  |CNzN(s)||2)ds 

where  we  have  used  the  Young’s  inequality.  From  (2.2),  we  have 
f”  M2 

||zN(t)!|2dt  <  -f  (w,  +  20  sup(|jCN||2  +  ||flN||2»||z||2 

for  all  z  e  Z.  Therefore,  the  theorem  follows  from  the  Datko’s  theorem  [7], 

(Q.E.D.) 

The  following  is  a  consequence  of  (3,  Theorem  6.9]  and  [2,  Theorem 

2.2]. 

Corollary  2.2  Suppose  ( A,B )  is  stabilizable  and  ( A,C )  is  detectable  and 
assume  (HI)  -  (H3)  are  satisfied.  Then  the  unique  nonnegative  solution  nN 
to  (1.8)  converges  strongly  to  n,  the  unique  solution  to  (1.5). 

Theorem  2.3  Suppose  that  B  is  compact  and  =  PNB  and  that 

(Hl)(i)  and  (H3)(i)  are  satisfied.  Then  (A,B)  is  stabilizable. 

Proof:  Let  us  consider  the  case  C  =  I  and  CN  =  PN  with  Y  =  Z.  Then  it  is 
easy  to  show  that  ( A.C )  is  detectable  and  (AN,CN),  n  *  1  are  uniformly 
detectable  since  (Hl)(i)  implies  that  for  some  M  *  1  and  w  independent  of  N, 


||SN(t)PN||  <  Me1**,  t  >  0.  It  then  follows  from  Theorem  1.1  and  (H3)(i)  that  for 
each  N,  (1.8)  with  CN  =  PN  has  a  unique  solution  nN.  Using  the  same 
argument  as  in  the  proof  of  Theorem  2.1,  we  have  ||nN||  <  S  for  some  positive 

t 

|  constant  0.  Thus,  by  Theorem  6.5  in  [3],  there  exists  a  subsequence  of  nN 

l 

!  converges  weakly  to  some  nonnegative,  self-adjoint  operator  ri  We  will  show 

|  that  n  satisfies  (1.5)  with  C  =  I.  Since  PNnN  =  nN,  BN*  riN 

0*PNnN  =  B*ff*.  Since  B*  is  compact,  for  each  z  €  Z,  Z?Nj  n^PNz  converges 

A  A 

j  strongly  to  B*Uz.  It  now  follows  from  [3,  Theorem  6.7]  that  n  satisfies  (1.5) 

,  But  since  (A,C)  is  detectable,  by  [10,  Theorem  3.2],  A  -  BB*h  generates  a 

uniformly  exponentially  stable  semigroup  on  Z. 

(Q.E.D.) 

Remark  2.4  Roughly  speaking.  Theorem  2.3  means  that  the  uniform 

stabilizability  implies  the  stabilizability  of  ( A,B ).  The  dual  statement  of 
Theorem  2.3  also  holds:  i.e.,  suppose  C  is  compact,  CN  =  CPN,  then  (Hl)(ii) 
and  (H3)(ii)  imply  that  ( A,C )  is  detectable.  This  statement  can  be  proved  by 
applying  the  exactly  same  arguments  as  in  the  proof  of  Theorem  2.3  to  the 
dual  Riccati  equation 

(AE  +  E A*  -  EC*CE  +  I)  z  =  0 


for  all  z  c  dom(/l*). 
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3.  Convergence  Rate 

In  this  section,  we  assume  that  (HI)  and  (H3)  are  satisfied  and  let 
=  pNs  and  CN  =  CPN.  Moreover,  we  assume 
(H4)  For  each  z  6  Z,  IIz  €  dom(/4*)  and  B  is  compact. 

From  (1.5),  we  have  for  all  z  e  dom (A) 

2<riz,Tz>  -  <B*m,B*nz>  +  <Cz,Cz>  =  0 
Thus,  (H4)  and  the  density  of  dom(/I)  in  Z  imply  that  for  all  z  £  Z 

2<A*nz,z>  -  <B*nz,5*nz>  +  <Cz,Cz>  =  0  . 

Define  the  self-adjoint  operator  nN  =  PNnPN.  Then  for  all  x  €  ZN 

(3.1)  2<rN*nNx,x>  -  <s*riNx,  5*rlNx>  +  <cNx,cNx>  +  <anx,x>  =  o  , 
where  AN  e  f(ZN)  is  a  self-adjoint  operator  defined  by 

(3.2)  <Anx,x>  =  2<(A*  -  AN*PN)Tlx,x> 

+  <B*( nN  -  nN)x,fi*(nN  +  nN)x>  for  all  X  e  ZN  . 

From  (1.8),  for  x  €  ZN 

(3.3)  2</tN*nNx,x>  -  <BN*nNx,£N*nNx>  +  <CNx,CNx>  =  0  . 

Hence  by  subtracting  (3.1)  from  (3.3) 

2<(tn  -  flNffN*nN)x,  (nN  -  rrN)x> 

+  <sN*(nN  -  nN)x,  fiN*(nN  -  nN)x  >  -  <anx,x>  =  o 

for  all  x  £  ZN.  Or  equivalently 

(3.4)  nN  -  nN  =  j"  e(JN-BNBN*nN)t((ItN  .  nN)sNsN*(nN  .  tf)  .  AN) 

x  °MN-*N*NV)t  dt 


I 


Similarly,  subtracting  (3.3)  from  (3.1),  wc  obtain 


(3.5)  nN  -  nN  =  f  c(^N-BN*N*nN)t((nN  .  nN)sNflN*(nN  .  ^  +  an} 

»  n 


x  e(/iN-BNBN*nN)tdt  . 


Here,  from  Theorem  2.1,  we  have 

||e^  B  n  PN||  S  M3e  3  .  t  5  0 


with  M3  ?  1  and  w3  >  0.  Since  nNPN  -  n,  strongly  by  Corollary  2.2  and 
B  is  compact 


||SN*(nN  -  nN)||  =  j|(nN  -  n)PNfl||  -  0  as  N  -  »  . 

Hence  by  the  variation  of  constants  formula  and  the  Gronwall’s  lemma, 

l|0U"-sNe',*n">t  pN||  ,  MsC(-“3*||sl|llsN*(nN-"N)||» 


It  then  follows  that  there  exists  an  integer  N0  such  that  if  N  J  NQ, 

* 

||e(y*  ~B  B  n  PN||  <  Mse  *  ,  t  *  0  . 


Now, 

from  (3.4) 

for  all 

X  € 

ZN 

(3.6) 

<(nN  - 

nN)x,x  > 

-r 

j(5N(nN 

■  AN)eMN-fiNfiN*nN)tx!|2  dt 

■-J 

.09 

<e(^N 

-fiNsN* nN)tXi  an  e(A*-BNBN* 

n%>dt 

and  from  (3.5) 

0 

(3.7) 

<(nN  - 

nN)x,x> 

-c 

|jfiN(nN 

-  ftN)c(x4N-BNBN*fiN)tx||idt 

-r 

Jo 

<e^N'» 

NBN*nN) tx  an e(AN-BNBN*nN 

>lx>dt  . 

These  inequalities  imply  that  for  x  t  ZN 


J 


- .  »'  '  ir  iV  !,»•/*  . .;*>*«'*•**  V** S  v&\V  ,V- '  ^M,v*-V,V.*  •*•****%’  *.  4*.* &,% / .****»  •,- V^V*.  .»/%.** ,4*v ’4y.*vi,. »f-‘ 
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Examples 

In  this  section  we  discuss  the  examples  in  which  (H1)~(H4)  are 


satisfied  and  then  apply  the  formula  (3.8)  and  (3.9)  to  obtain  a  convergence 
rate  of  nN  to  IL 

4.1  Parabolic  Systems 

Assume  V  and  H  are  Hilbert  spaces  and  V  c  H  with 

continuous  dense  injection  i.  Consider  a  bilinear  form  o  on  V  such  that 

(4.1)  |o(u,v)|  S  c||u||  ||v||  for  u,v  €  V 

v  v 

(4.2)  o(u,u)  ?  w|ju||^-  P||u||2  for  u  «  V 

where  u  >  0.  It  then  follows  from  [9]  that  there  exists  an  operator  A  € 
*(V,V*)  such  that 

(4.3)  o(u,v)  =  <-^u,v>v+v  for  u,v  €  V  , 

where  V  C  H  =  H*  C  V*  and  H  being  the  pivoting  space,  and  that  A 
on  H  with 

(4.4)  dom(y4)  =  (x  e  H  :  Ax  e  H)  dense  in  V, 

generates  the  analytic  semigroup  on  H  and  V*.  For  given  B  e  £(U,H)  and 
C  e  X(H,V)  consider  approximating  problems  (ZN,/iN,SN,CN);  i.e.  let  ZN  be  a 
sequence  of  finite  dimensional  subspace  of  V  and  ^N:  ZN  -*  ZN  is  defined  by 

(4.5)  <-/lNz,x>  =  o(z,x)  for  z,x  e  ZN  . 

Let  PN  be  the  orthogonal  projection  of  H  onto  ZN  and  assume  BN  = 
PnB  and  CN  =  CPN.  We  assume  the  approximation  condition: 
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For  each  z  €  V  there  exists  an  element  zN  €  ZN 


(4.6)  such  that  jjz  -  zN||  <  e(N)  where  «(N)  -*  0  as  N  -  “  . 

v 


It  then  follows  from  [2]  that  (HI)  holds  and  if  (A.B)  is  stabilizable  and  (A,C)  is 
detectable,  then  (H3)  holds.  Thus  from  Corollary  2.2,  nN  converges  strongly  to 
II  However  one  cannot  apply  the  formula  (3.8)-(3.9)  as  it  is,  since 
nZ  C  dom(,4*)  is  the  maximal  regularity  without  assuming  any  regularity  of  C. 
This  can  be  demonstrated  by  the  following  example.  Consider  the  case  when 
H  =  L2(0,1)  and  V  =  H*(0,1),  and 


f1  d  d  , 

o(u,v)  =  -  u(x)  - v(x)dx  for  u,v  e  Hi  . 

Jn  dx  dx  u 


Let  us  consider  the  Liapunov  equation  on  H 


(4.7)  AE  +  EA  +  Q  =  0 

where  Q  is  self-adjoint  operator  on  H.  If  for  each  z  e  Z, 


(Ez)(x)  =  f  <Kx,y)z(y)dy  and  (Qz)(x)  =  1  q(x,y)z(y)dy, 
•  0  ^  0 


=  [  q(x,y 

J  n 


then 

<p  satisfies  A0  h 

J-  q 

=  0  with 

Dirichlet  boundary 

a2 

a2 

where 

+ 

-©■ 

n 

< 

<p  for  Q 

€  H2([0,I]  x  [0,1]). 

3x2 

ay2 

[6], [8]) 

f  f  f|-* 

Jo  Jo  J'Sx1 

2 

+ 

i -*n 

3y 2  1  J 

dxdy  S  M  f  f  |q| 
•*o  ''O 

dxdy  . 


This  implies  EL2  C  dom(-4)  is  the  maximal  regularity. 

Hence  we  will  modify  the  arguments  in  Section  3  to  improve  the 
formula  (3.8)-(3.9)  for  this  example.  First  we  note  that  in  (3.2)  for  x  e  ZN 


| <(/i*-^N*pN) nx,x>j  =  |c*x,(nN  -  n)x)| 

<  cIIMIJI<n  •  ""I'llv  **  <4» 


C  •*  .  '•  ’  4  \  »  5 .  *  *  ?  *  ‘  •’  -  *  "  •-  *  A*  ■*»  •  '!  S4I  •  **  .**  *S  -  *:  .  ■*  .  /  '  ‘ 


Thus  from  (3.2) 


(4.8,  7  -  sup  ig^L  (2c  sup  ll<n-fiK)X  llv 


I*  II  *>  Ml 

\r  » 


x*0 


[|g*(n-nN)x  [|v 


+  2a0|jfl||  sup 

x*0  llx 


where  ||x||h  <  «||x||v. 


Lemma  4.1 


There  exists  a  positive  constant  M  such  that 


f°°  liN  DNRN*nNw  M  2  2 

jje^  8  n  PNx|j  dt  <  Mjjx ||  ,  and 

J0  V  H 

I-”  /  jN  RNRN*nNNt  M  2  2 

J  B'(  ’  p  >llv  dI  £  Mllxl!„  ' 


Proof:  Let  (N(t)  =  e^x  B  n  )t  pNx  ,  >  Q  Then  fN(t)  satisfies 


dt 


$N(t)  =  (,4N-flNflN*nNKN(t)  ,  t  >  0  , 


so  that  from  (4.5) 


t  l'H(t)£  +  °<^N)  = 


</(N(t),BN*nN(N(t)> 


H 


and  from  (4.2) 


1  d 


y  d r  iisN(°Hh  *  "ll'N<,)llv  5  <p  +  2Bllsll >  ll!N(,)ll„  • 

The  integration  of  this  inequality  with  respect  to  t  yields 

y  II'n«>|4*  “  fQ  PN<»C  <  y  ll'N<°>l|I  ♦  <f  ♦  2Bllsll’>  fo  P”<s)|Cd!- 

Now  the  lemma  follows  from  Theorem  2.1. 


Q.E.D. 
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It  then  follows  from  Lemma  4.1  and  (4.8)  that 
f  <an  e(^N-5N5N*nN)tx  f  e(^N-BNsN*nN)tx  > 


<  7  |  ||e^N'^N^N  nN)f  P*x„*  dt  <  M7||x||^  . 

/  iN  dNdN^Nv, 

Similary  for  e^  '  ,  t  >  0.  Therefore  we  obtain,  using  (3.6)  and 


(3.7), 

(4.9)  ||nN  -  nN||  <  M7  . 
where  7  is  given  by  (4.8). 

Consider  the  (1-dimensional)  parabolic  control  system  [2); 


d  3d  0 

—  z(t,x)  =  —  (p(x)  —  z)  +  q(x)  —  z  +  r(x)z 

ot  3x  ox  dx 

m 

+  E  bj(x)u;(t)  in  (0.1) 

i=i 

with  boundary  condition  z(t,0)  «  z(t,l)  =  0,  where  p  £  C^O.l),  being  bounded 

d  „ 

below  by  a  positive  constant  w,  —  q,  r  €  L  (0,1),  and  b:  £  L2(0,1),  i  =  1 . m. 

dx  i 

In  this  case,  H  =  L2(0,1)  and  V  *=  Hq(0,1),  and  the  bilinear  form  a  is 
given  by 


o(u,v) 


d  d  d 

—  u  —v  -  (q(x)  _u  +  r(x)u)v)dx  . 
dx  dx  dx 


B:  R"1  -*  L2(0,l)  is  defined  by 


(Bu)(x)  =  E  bj(x)uj  for  u  €  R"1  , 
i=l 

and  dom(/l)  =  dom(^*)  =  H2(0,1)  O  Hq(0,1).  Let  us  consider  the  following 
finite  dimensional  subspace  ZN  of  V: 

N-l 

ZN  =  {z  £  H  :  z(x)  =  E  a;Bf(x)  ,  cc;  £  R) 

i=l 


15- 


r 


where  B?V).  i  =  1.....N-1  are  the  linear  B-spline  elements  on  the  interval  [0,1]; 


i.e.. 


i+1 

-N(x- - )  , 

N 

N(x- -Li- )  , 
N 


x  e 


x  e 


&-t] 

K-lr] 


0  ,  elsewhere 

Then  the  approximation  condition  (4.6)  is  satisfied  [8].  Suppose  (A,B)  is 
stabilizable  and  (A,C)  is  detectable.  Then  (1.5)  has  the  unique  solution  n 
and  using  a  similar  arguments  to  those  given  above  to  show  the  regularity  of 
solutions  to  Liapanov  equation  (4.7),  one  can  show  that  for  x  e  H,  Fix  e 
dom(/l*).  Since  A*  is  closed  in  H  and  dom(/i*)  C  V,  by  the  closed  graph 
theorem,  there  exists  a  positive  constant  kj,  such  that  ||nz ||H2(0  jj  4  killzllL2(0i) ' 
Hence  the  fundamental  error  estimate  (e.g.,  [8])  gives 

llnz  •  "NzIIli  <  k>  (r)!HzIIli 
-  "Nzll„;  ‘  k»  [f;]llzllLi 

for  some  positive  constants  k2,k3.  Now  it  follows  from  (4.8)  and  (4.9)  that 
|jnN  -  nN ||  <  k  j  for  some  constant  k. 

4.2  Hereditary  Differential  Systems 

Consider  the  hereditary  differential  system  in 

x(t)  =  A0x(t)  +  AjX(t-r)  +  P  A(8)x(t+8)d8  +  Bu(t) 

J  -r 

(4.10)  {  x(0)  =  n  and  x(9)  -  #8),  -r  <  8  <  0 

y(t)  =  C  x(t) 


-16- 


•ya 


and  the  optimal  control  problem;  for  given  initial  data  z  =  (b,<>)  £  R"  x 


L2(-r,0;  R")  ,  minimize  the  cost  functional 


(4.11)  J(u,z)  =  j  (|y(t)|l+  ju(t)jJ)dt  . 

*  r\ 


Here,  x  £  R",  u  e  R"1  and  y  t  R*1  and  the  element  of  A()  is  square 


integrable.  It  is  shown  [1]  that  (4.10)  and  (4.11)  are  equivalently  formulated  as 


the  problem  (1.1)  -  (1.4)  on  the  product  space  Z  =  R"  x  L2(-r,0;R");  i.e.,  z(t)  = 


(x(t),x(t+ •))  e  Z  is  the  mild  solution  of  (1.1)  with 


domU)  =  {(b,<#>)  €  Z  :  <f>  £  H^-lO)  and  <K0)  =  b}  , 


for  (<K0),<f>)  £  dom(/l) 


=  (Ao0(O)  +  A,0(-r)  +  A(8)0(0)d9,  0) 


The  input  operator  B  :  R"1  -*  Z  and  the  output  operator  C  :  Z  -*  R1*  are 


given  by 


Bu  =  (Bu,0)  e  Z  and  C(n,<J>)  =  Cb  . 


Let  us  consider  the  averaging  approximation  [1]  of  (4.10);  let 


ZN  =  (z  £  Z  :z=  (a0,  Y.  akx  ;  j.,  ),  ak  €  Rn  ,  0  <  k  <  N)  C  Z  , 


k=l  [-— r>-  —  r) 
N  N 


and  /1N  has  the  matrix  representation  (Qn)_1Hn  on  Rn(N+0  when  ZN  is 
identified  with  Rn(N+|)  by  its  coordinate  vector  col(aJ,...,aJ),  where  the  block 
diagonal  matrix  QN  and  the  block  Hessenberg  matrix  HN  are  given  by 


rAN  an  an 


qn  =  N 


and  HN  = 


I  -I 


I  -I 


-17- 


wit  h  Aq  =  AQ 


i-1  N-l 

,  AjN=  |  nA  (0)d0  and  Ajj  =  Ax  +  |  N  A(0)d0.  Note  that 


PNfi  =  B  and  CPN  =  C.  Set  BN  =  B  and  CN  =  C.  Then  (/4N)*  has  the 
matrix  representation  (QN)_1HN  on  pniN+i),  (Hl)(i)  is  proved  in  [1]  and 

(HI)(ii)  is  proved  in  [3].  Using  the  arguments  in  [5],  [7]  one  can  show  that 
(H3)  is  satisfied  (i.e.,  (i)  is  straightforward  but  (ii)  is  not  so).  Thus,  the 
formula  (3.8)-(3.9)  yields 


||nN  -  nN||  «  2 1| (A*  -  /iN*PN)n|| 

By  the  regularity  result  in  [4],  if  A(  )  e  H^-r.OjR"*"),  then 
A* n  +  C*C  «  dom(,0 

where  dom(/l*)  =  {(y,0)  €  Z  :  0  e  H1  and  0(-r)  =  Ajy}  and  A*( y,0)  = 

(0(0)  +  ^Jy,  -  m  +  ^T(-)y)  «  Z  [3],  Since  C*C(n,0)  =  (CTCn,0)  t  Z  for 

(n,0)  €  Z,  this  implies  that  if  II z  =  (y ,0),  then  0  €  H1  so  that  0  e  H2,  and 

since  A*  is  closed,  ||0||  <  M||z||  f°r  some  constant  M.  It  then  follows 

H2  z 

from  the  arguments  and  error  estimate  in  [1],  [3]  that 
„  M 

||(-4*-^N  PN)n(y,0)||  <  — (|y |  +  ||0 1|  2) 

✓N  h 

Hence  we  obtain  ||nN  -  nN||  =  0(  — )  . 
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